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Abstract 

Knowledge about the location of licensed primary-users (PU) could enable several key features in 
cognitive radio (CR) networks including improved spatio-temporal sensing, intelligent location-aware 
routing, as well as aiding spectrum policy enforcement. In this paper we consider the achievable accuracy 
of PU localization algorithms that jointly utilize received-signal-strength (RSS) and direction-of-arrival 
(DoA) measurements by evaluating the Cramer-Rao Bound (CRB). Previous works evaluate the CRB for 
RSS-only and DoA-only localization algorithms separately and assume DoA estimation error variance 
is a fixed constant or rather independent of RSS. We derive the CRB for joint RSS/DoA-based PU 
localization algorithms based on the mathematical model of DoA estimation error variance as a function 
of RSS, for a given CR placement. The bound is compared with practical localization algorithms and the 
impact of several key parameters, such as number of nodes, number of antennas and samples, channel 
shadowing variance and correlation distance, on the achievable accuracy are thoroughly analyzed and 
discussed. We also derive the closed-form asymptotic CRB for uniform random CR placement, and 
perform theoretical and numerical studies on the required number of CRs such that the asymptotic CRB 
tightly approximates the numerical integration of the CRB for a given placement. 
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Cramer-Rao Bounds for Joint RSS/DoA-Based 
Primary-User Localization in Cognitive Radio 

Networks 

I. Introduction 

Cognitive Radio (CR) is a promising approach to efficiently utilize the scarce RF spectrum 
resources Q2. In this paradigm, knowledge about spectrum occupancy in time, frequency, and 
space that is both accurate and timely is crucial in allowing CR networks to opportunistically 
use the spectrum and avoid interference to a primary user (PU) [0. In particular, information 
about PU location could enable several key capabilities in CR networks including improved 
spatio-temporal sensing, intelligent location-aware routing, as well as aiding spectrum policy 
enforcement 0. 

The PU localization problem in CR networks is in general different from localization in other 
applications such as Wireless Sensor Networks (WSN) flU and Global Positioning System (GPS) 
0, due to the following two features. First, a PU does not communicate with CRs and very 
limited knowledge about its signaling, such as transmit power or modulation scheme, is available. 
Therefore, passive localization techniques should be applied. Second, since CRs need to detect 
and localize PUs in the whole coverage area at a very low signal-to-noise ratio (SNR), in order 
to avoid interference to the primary network, the required number of CRs is relatively large and 
cooperation among CRs is necessary. 

A. Related Work 

Prior research on passive localization can be categorized into three classes based on the 
types of measurements shared among sensors to obtain location estimates [4]. Received- signal- 
strength (RSS) based algorithms use measured received power from the PU to provide coarse- 
grained estimates at a low hardware and computational cost. Time-difference-of-arrival (TDoA) 
based algorithms obtain location estimates from time differences among multiple receptions of 
the transmitted signal. They are not suitable for CR applications since TDoA-based algorithms 
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require perfect synchronization among CRs. Direction-of-arrival (Do A) based algorithms use 
target DoA estimates observed at different receivers to obtain location estimates. The DoA 
estimates can be obtained from either multiple antenna arrays, directional antennas or virtual 
arrays formed by cooperative CRs. Therefore, RSS and DoA-based algorithms are the proper 
choices for the PU localization problem. 

Passive localization based on RSS and/or DoA measurements has been well-studied in both 
WSN and CR literature, in terms of algorithm design and evaluation. RSS-based algorithms 
include range-based and range-free techniques, a detailed survey can be found in O. Popular 
DoA fusion algorithms include Maximum Likelihood (ML) 0, Stansfield algorithm (HI and 
Linear Least Square (LLS), each of which provides different tradeoff between accuracy and 
complexity. Algorithms that jointly utilize RSS and DoA information are also proposed. For 
example, weighted versions of the popular DoA fusion algorithms are shown to improve the 
localization accuracy J7J. The weights are typically estimated error variances of individual DoA 
measurements, that are obtained from instant RSS |[9]|. 

In this paper, we focus on characterizing the achievable performance of joint RSS/DoA-based 
localization algorithms by means of the Cramer-Rao Bound (CRB), which provides a lower 
bound on the estimation accuracy of any unbiased estimator using either or both RSS and 
DoA. The CRB for RSS-only localization is studied in [31 and [[Toll assuming independent and 
correlated shadowing channel, respectively. The DoA-only CRB is discussed in several papers 
0, lHTfl - lfl4l . however they all assume the DoA estimates from different CRs are subject to 
i.i.d. Gaussian errors with zero mean and fixed variance. Since the DoA estimation error variance 
depends on instant RSS and other system parameters (e.g. number of antennas and samples, array 
orientation error) [9], assuming it's constant makes the derived bounds less accurate. Recently, 
the mathematical model of DoA estimation error is used in lfT5l - lfT71 to calculate DoA estimation 
error variance for some specific system settings, such as node density and placement, using fixed 
value of RSS for all CRs. The variance values are then plugged into the existing DoA-only CRB 
to evaluate the localization performance. This approach does not include the system parameters 
in the mathematical model of DoA error variance and channel parameters in RSS, thus the CRB 
is less comprehensive. 
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B. Contributions 

Since the DoA estimation accuracy strongly depends on RSS, and practical algorithms show 
that using both RSS and DoA estimates provides better localization accuracy than using DoA 
alone, the CRB for joint RSS/DoA-based localization is more capable to characterize the posi- 
tioning ability of the CR network. Along this direction, this paper has the following contributions: 

1) Derivation of the CRB for joint RSS/DoA-based passive localization considering the 
interdependency between RSS and DoA, for a given CR placement. The CRB is derived for 
both optimal DoA estimation algorithm, with error given by the CRB of DoA estimation 
error variance, and MUSIC algorithm. Therefore, the derived CRB provides both ultimate 
achievable accuracy and the accuracy achieved by a practical DoA estimation algorithm. 

2) Derivation of closed-form asymptotic joint RSS/DoA-based CRB for uniform CR place- 
ment with i.i.d shadowing. The asymptotic RSS-only CRB for such case is also derived 
as a byproduct. We also provide theoretical and numerical studies on the required number 
of CRs that ensures the asymptotic CRB tightly approaches the numerical integration of 
the CRB for a given placement. 

3) Study of the impact of various system, channel and array parameters on the localization 
accuracy, theoretically and via simulation. The considered parameters include number of 
CRs, shadowing variance, correlation distance, number of antennas, number of samples 
and array orientation error. The study provides insights and guidelines for practical system 
design of PU localization in CR network. 

4) Comparison of the derived CRBs with practical RSS-only and joint RSS/DoA localization 
algorithms. We study the robustness of practical algorithms to major system parameters, 
and reveal the potential improvements of practical localization system and algorithm design. 

The rest of the paper is organized as follows. The system model is introduced in Section|IIl The 
joint CRB derivation for a fixed CR placement is presented in Section Unl The asymptotic joint 
CRB derivation for uniform CR placement is shown in Section [IVJ Numerical results evaluating 
impacts of various parameters on the joint CRB are discussed in Section |V] Finally, the paper 
is concluded in Section EO 
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II. System Model 

Assume N CRs cooperate to localize a single PU. 2-dimensional locations of the PU and the 
n th CR are denoted as £p = [xp,yp] T and £ n = [x n ,y n ] T , respectively. All locations are static 
in the observation period, and CR locations are known. Available measurements at CRs are RSS 
and DoA. RSS measurement at the n th CR is modeled as ib n = P T coW ^" 7 10 Watt, where P T 
is the PU transmit power, c is the (constant) average multiplicative gain at reference distance, 
d n — \\£ n — £p\\ is the distance between the n th CR and PU, 7 is the path loss exponent, and 
10 - W 10 i s a random variable that reflects shadowing. The RSS is normally expressed in dBm 
using the transformation <f> n = 10 log 10 (100Ch/> n ), the result is given as 

4> n = 10 log 10 (1000P T c ) - IO7 log 10 d n -s n = ~4> n - s n , dBm. (1) 

Denote collection of RSS measurements from all CRs as 4> — [4>u 02, • • • , 4>n] t - The conditional 
distribution of 4> (f° r a given £ P ) is <fi ~ A/"(0, fl s ), where = [<f) 1 , <p 2 , . . . , 4> N ] T , and f2 s is 
the covariance matrix of the collections of the shadowing variables s = [s±, S2, ■ ■ ■ , sn] t given 
by {fl s } mn = afe-Wtm-inW/Xc^ w j iere x c is the correlation distance within which the shadowing 
effects among nodes are correlated. 

The DoA of the PU at the n th CR is given as 9 n = arctan(^5^-) = Z(£ P ,£ n ). CRs 
perform array signal processing techniques, such as MUSIC lfT8l or ESPRIT |fT9l , to obtain DoA 
estimates. The estimated DoA is commonly modeled as 9 n = 6 n + v n flU, where v n ~ A/"(0, a^) 
and is the DoA estimation error variance. We denote 9 = [81, 82, ■ ■ ■ , 8n] as collections of 
DoA measurements of all CRs at the fusion center. We consider two different modelings of the 
DoA estimation error variance, using CRB and MUSIC algorithm, respectively. The CRB of 
DoA estimation error variance for unbiased DoA estimators using Uniform Linear Array (ULA) 
is given by 10 

n ' CRB ~ (Kcos8 n yN s N a (N*-lW ( ) 

where n is a constant determined by the signal wavelength and array spacing, N s is the number 
of samples, N a is the number of antennas, 8 n is the array orientation with respect to the incoming 
DoA defined as 8 n = 8 n — 8 n , where 8 n is the orientation of the n th ULA, and p n is the signal- 
to-noise (SNR) ratio given by p n = ip n /PM, where Pm is the measurement noise power at the 



5 



71 th CR. Using the definition of SNR we can simplify © as 

2 6Pm 1 1 _ at ,2 \ 1 /on 

where f3 = _ 2M tfV M 2_u an ^ fcRB(4>n) — t - . The DoA estimates from MUSIC are, asymptot- 
ically in the sample size, unbiased and Gaussian distributed JH. The estimation error variance 
using ULA is given by O 

16 1 1 

(k cos # n ) 2 N s N a (N£ - l)p n N a p n cos 2 Q n 

where fuuQn) - '^ n+i f 2 l/Na) . Note that both © and © depend on RSS and PU location (used 
in calculation of 9 n ). 

We do not assume any specific distribution of CR placements in the derivations of Sec. HH] 
For the asymptotic CRB derived in Sec. [IV] we assume the CRs that can hear the PU form a 
circle with radius R and are uniformly placed in the area. For this scenario, the distribution of 

9 n is given by 9 n ~ U[0, 2ir), and the distribution of d n is derived as 

2r 

Ro<r<R (5a) 



p dn {r)= { (R 2 -R 2 oV 

0, otherwise, (5b) 

where R is a guard distance to avoid overlap between CRs and PU. The CRs are placed 

independently within the area, which indicates independency among all re 's, all cZ n 's, and pairs 

of 9 n and d n . We further assume the array orientation error is distributed as 9 n = 9 n — 9 n ~ 

U{—9 T ,9 T ), which means the array orientation is uniformly distributed around the incoming 

DoA with parameter 9 T . The PU localization problem uses RSS and DoA measurements to 

obtain PU location estimate £p = [x p , y p ] T . A brief table of notations used in the paper is given 

in Table E 

III. Joint CRBs for Fixed CR Placement 

In this section we derive the joint CRB and the corresponding bound on RMSE for a fixed 
CR placement, for DoA estimates obtained from both optimal estimator and MUSIC algorithm. 
We also derive the CRB and RMSE for RSS-only localization as a byproduct. Using RSS and 
DoA as measurements, the covariance matrix of unbiased estimation of PU locations £p is 
lower-bounded by the CRB 
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where F is the 2x2 Fisher Information Matrix (FIM) given by 

d 2 



E 



e,<p 



d£ 2 P 



\ogp(e,c/>\ep) 



(7) 



Therefore the RMSE is bounded by RMSE > a/{F _1 }h + {F _1 } 2 2, where {X}^. denotes the 
ijth e i emen t of matrx X. Using the standard decomposition of conditional probability p(6, 4>\£p) = 

p(6\(f),£p)p(4>\£p), the FIM is decomposed as 

d 2 

— T \ogp(0\(j),£ P 



E 



e,4> 



d£i 



E- 



d 2 



d£ 2 P 



logp{<f>\£i 



4>- 



(8) 



Note that F^ is the FIM for using only RSS to localize the PU, therefore its inverse bounds the 
localization accuracy of algorithms that only use RSS readings. In the rest of the section, we 
first derive the RSS-only FIM F^. We then present the results for the joint FIM F by deriving 



e\<t> 



for optimal DoA estimator and MUSIC algorithm, respectively. 



A. RSS-only CRB 

To derive the RSS-only FIM F-?, we first explicitly express the logarithm of the PDF of 4> 



\ogp(ct>\£p) = 
The RSS-only FIM F^ is then given by 



log [(2vr) Af / 2 (deta) 1 / 2 ] - l -Q>- c/>fn;\4> - 4>). 



2 * 



dtp 



(9) 



(10) 



The elements of F-r are derived as 

4> 



ii 



22 



1.2 



d{<f> - W d{ct> - <t>) = 2 Ax r D - 2f2 i D - 2Ax 

OXp OXp 



-n: 



4> 



dy P dy P 

21 



e 7 2 Ay T D- 2 n g - 1 D- 2 Ay 



e 7 2 Ax r D- 2 a- 1 D- 2 Ay, 



(11) 

d N ), 



dxp dyp 

where e = 100/(log 10) 2 , and vectors and matrices are defined as D = diag(<ii, <i 2 , 
Ax = [Azi, Ax 2 , . . . , Ax N ] T , Ay = [Ay u Ay 2 , . . . , Ay N } T , and Ax n = x P - x n and Ay n = 
Up — Vn- Detailed derivations are provided in Appendix A. To obtain a compact expression of 
Fg, let's define L = [Ax, Ay] T and A = ^2-D 2 fi s D 2 . Then it is straightforward to verify that 



LA L T . Therefore, the RMSE of RSS-only PU localization is bounded by 



RMSE R , F > ^{FZ 1 }!! + {F^} 22 , 
where the subscript R, F stands for RSS-only bound for fixed placement. 



(12) 
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B. Joint CRB using Optimal DoA Estimator 

In this subsection we derive the joint CRB with DoA estimations given by the optimal 
estimator, using the DoA error variance given by <r 2 CRB . To derive the conditional FIM of 
DoA given RSS Fg^, we first explicitly express logarithm of the conditional PDF p(9\4>,£p) as 



N 



\ogp(e\ct>,£ P : 

Then Fg^ is given by 



^2 \ log(cos 9 n ) - ^ lc 

n=l 



27T/3/, 



CRBVP n I 



COS 2 9 r 



n ) 2 



2/3/, 



(13) 



crbW 



N 



0\4> 



EN 



where g n 



n=l 

A 



d 2 9n 
d£ 2 P 



E 



e,cp 



d 2 K 
di 2 P 



(14) 



= c °l^ n B ^ and h n = log(cos^). The elements of Fg^ are derived as 



+ 2 tan 2 9 n 



9\4> 



o\<t> 



22 
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-t¥\ 

n=l a ™ 


a cos 2 9 n 


dl 


n=l fl ™ 


a cos 2 9 n 


[ dl 


= { F§ l^21 


n=l 



+ 2 tan 2 9 r 



Ax n Ay n acos 2 9 n 



dl 



dl 



2 tan 2 9 n 



(15) 



where a = cqPtg^^ 2 ^ / (3. Detailed derivations are provided in Appendix A. To obtain a compact 
expression of Fg,^, let's define P = [Ay, — Ax] T and T = diag(7i,7 2 , . . . ,Jn), where 7„ = 
j_ | °cm 2 o„ + 2tan 2 n j. Then it is straightforward to verify that Fg^ = PrP T Therefore, the 
joint FIM and the corresponding RMSE are given by V JF c = PTP T + LA _1 L T and 



RMSE J>FiC > J{Vj* F>c } n + {FJ>, C } 



22- 



(16) 



where the subscript J, F, C represents joint CRB for fixed placement using CRB of DoA 
estimation error variance. 



C. Joint CRB using MUSIC Algorithm 

In this section we derive the joint CRB with DoA estimations given by the MUSIC algorithm, 
using the error variance given by cr 2 MU . The conditional FIM of DoA given RSS Fg.j is derived 
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by replacing fcRB{4>n) in (TT3T ) with fMui^n) given by ©. Applying the results in Appendix A, 



we obtain Fg^ = PAP , where A = diag(<5i, 6 2 , ■ ■ ■ , 5 N ) and 



cos 2 9 % 



CoP^/M _ Pm 



P 



11 



E- 




+ 2 tan 2 n 



(17) 



Since the RSS-only FIM is independent of DoA estimation algorithm, the joint FIM and the 
corresponding RMSE using MUSIC algorithm are given by Fj^m = PAP T + LA _1 L T and 

RMSE J)FjM > V /{F7> iA/ } 11 + {F7> iM } 22 , (18) 

where the subscript J, F, M represents joint CRB for fixed placement using MUSIC error 
variance. Note that the above calculations are based on a particular fixed node placement. If 
the nodes are randomly placed within the area according to some distribution, say uniform 
distribution, integrating (fl"2~l) . (fT6l) and (fl"8l for 2N times over all node coordinates gives the 
theoretical average performance bound. 

IV. Asymptotic Joint CRB for Uniform Random CR Placement 

The bounds on RMSE derived in Section [HI] are useful to evaluate the achievable localization 
performance for a given CR placement. Numerical integration or ensemble averaging need to be 
performed to get the average accuracy for random CR placements. In this section we characterize 
the achievable localization accuracy for the most common random placement, i.e. the uniform 
random placement, by deriving the closed-form asymptotic CRB for such case, assuming i.i.d. 
shadowing and optimal DoA estimator. We also provide theoretical analysis on the required 
number of CRs as a function of system parameters to make the asymptotic bound tight. 

A. Asymptotic RSS-only CRB 

For the RSS-only FIM F^ with i.i.d. shadowing, we first rewrite it from CCD) as 



1f~ 

N + 



i2_ 



N 



Y, d n 



n=l 



COS 6 n 

sin^ 



Icos 9 n , sin 6 r , 



(19) 



in which we use the fact that [Ax n , Ay n ] = [d n cos n , d n sin 9 n ). From (fT9l) we observe that 
-^F^ can be interpreted as the ensemble average of a function of random variables d n and 9 n , 
with statistical mean given by 



1e 

N 



2at 



67 2 log(i?/i?o 

a*(R 2 ~ R 2 ) 



(20) 
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where we obtain the first equality from (fl9l) based on the i.i.d. distribution of d n and 9 n , and 

E 6n { [cos 9 n , sin 9 n ] T [cos9 n , sin# n ] } = |l 2 , where I 2 is a 2 x 2 identity matrix. The deviation 



probability of the ensemble average -^F^ from the statistical mean E 
following theorem. 

Theorem 1: (Deviation Probability of RSS-only FIM) For any 5 > 0, we have 



is given by the 



Pr 



If,-1e 



> 5 } < 



26 2 7 4 



2R 2 R 2 



lo^(R/R ) 
(R 2 



(21) 



where ||-|| represents Frobenius matrix norm. 

The proof of the theorem is based on Chebyshev's inequality for random matrix, and all details 
are provided in Appendix B. The intuition from Theorem Q] is that we can well-approximate ^F^ 



with jj E 



when N is large enough. In order for the deviation probability to be smaller than 
a pre-defined threshold rj E (0, 1), we need to bound the right-hand-side of (1211) by 77. The 
required number of CRs is then bounded by 



2„,4 



N > 



2R 2 R 2 



fo£(R/Ro 
(R 2 



Rl) 2 



(22) 



Applying the approximation, the RMSE of RSS-only algorithms for uniform random CR place- 
ment is given by 

2 



1/2 



RMSE^ > 



(23) 



where the subscript R, U represents RSS-only CRB for uniform placement. 



B. Asymptotic Joint CRB 

We apply similar procedure in Sec. IIV-AI to derive the asymptotic joint CRB from F 
F^ + Fg^. We first rewrite Fg^ from CES) as 



-F~r 

N e ^ N 



i N 

— [ ad n il+2) cos 2 9 n + 2d~ 2 tan 2 9 n 

71=1 



sin 9 n 
- cos 9 r 



[sin 9 n , — cos 9 j 



(24) 



From (l24l) we observe that ^Fg^ can be interpreted as the ensemble average of a function of 
random variables d n , 9 n and 9 n , with statistical mean given by 



N 



e\4> 



(R 2 - R 2 ) 



2hg(R/R ) 



'tan On 



a 



-(R-y-R^) 
27 



sin 20 5 
29 T 



+ 11 



2- 



= /g |$ (£,7,a s 2 ,0 T ,/?)I 2 . 



(25) 
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The deviation probability of the ensemble average -k F from the statistical mean jr E [F 



i E 

N ^ 



±- E 

N ^ 



0\4> 



is given by the following theorem. 



Theorem 2: (Deviation Probability of Joint FIM) For any <5 > 0, we have 



Pi 



N 



N 



E [Fl 



> So > < 



MS* 



6 2 7 4 



E[/„] 



2e 7 2 log(iVi? ) 
a*(R 2 ~ R 2 ) 



(26) 



where f n = a<in ( ' 7+2 ' ) cos 2 9 n 



2d~ 2 tan 2 n , E [f n ] = 2f d $(R, 7, a 2 , T , /3) and 



E [fn] 



4 tan 5 # T 



x 2 Fx 



5 -l 7 - 

2' '2' 



a 2 (sin 46> T + 8 sin 26> T 



- tan 2 9t 
369 T 



2«(^ (7+2) -^o (7+2) ) / sm2 T 



(7 + 2)(i? 2 



*8) 



7 T 



H-2(7+l) _ i?" 2 ^ 1 ) 



(27) 



32^ T ( 7 + 1)( J R 2 - J R 2 ) 
and 2-^1 (a;) is the hypergeometric function. 

The proof of the theorem is provided in Appendix B. The intuition from Theorem [2] is that 
we can well-approximate 4F with — E [F] when N is large enough. In order for the deviation 
probability to be smaller than rj, we need to bound the right-hand-side of (1261) by r). The required 
number of CRs is then bounded by 



N > 



e 2 7 4 

<R 2 Rl 



e [/a 



E[/„] 



2e7 2 log (i?//? ) 
^ 2 



n 2 



(28) 



Applying the approximation, the joint FIM and the corresponding RMSE are given by ~Ej } u,c 



N 



^( J R,7,a 2 ) + %( J R,7,a 2 ,^,/3) 



I 2 and 



1/2 



RMSE 



j,u,c 



> 



N 



(29) 



^( J R,7,a 2 ) + /g | ^( J R,7,a 2 ,^,/3) 

where the subscript J, U, C represents joint CRB for uniform placement using CRB of DoA 
estimation error variance. Note that different from the RMSE expression ([TBI which is condi- 
tioned on a fixed CR placement, the result in (1291 provides asymptotic closed-form expression of 
RMSE for uniform CR placement. Since the Chebyshev's inequality is not a very tight bound on 
the derivation probability, the results we obtained may overestimate the required number of CRs. 
Although similar analysis can be performed based on more accurate mathematical tools, such as 
central limit theorem (approximate the random matrix as Gaussian distributed) or Chernoff bound 
(requires the tail of the random variable distribution to be Gaussian-like) [20], the results obtained 



1 1 



by the Chebyshev's inequality are more general in that they do not assume any distribution of 
the random matrix. Therefore, (|22l) and (1281) mainly serve to provide insights of impact of system 
parameters on the required number of CRs to make the asymptotic bounds accurate. 

V. Numerical Results 

In this section we present numerical results to verify the derived CRBs and study the achievable 
localization performance for different node density, channel conditions and array parameters. We 
compare the derived bounds with the joint CRB with constant RSS derived in Eqn (6) of lfl6l to 
show the theoretical advantage of the new bound. Performance of practical localization algorithms 
is also included to investigate how close the CRBs can be achieved. A numerical study on the 
accuracy of the asymptotic CRBs derived in Section IV is presented at the end of the section, 
highlighting the conditions when the asymptotic CRBs tightly approximate the exact bounds. 
We start by a brief introduction of the considered localization algorithms. 



A. Practical Localization Algorithms 

We have selected two cooperative localization algorithms, one for RSS-only localization 
and the other one for joint RSS/DoA localization. For the RSS-only case, Weighted Centroid 
Localization (WCL) is a range-free algorithm that provides low-complexity and coarse-grained 
location estimates. In this technique, PU location is approximated as the weighted average of 
all secondary user positions within its transmission range, where the weights are proportional 
to RSS of each CR. The WCL location estimate is given by ipyvch = X^=i Vw^n/ S^Li ^n- 
WCL does not require knowledge of PU transmit power and channel conditions, and is robust 
to shadowing variance, compared to range-based algorithms such as Lateration 

For the joint RSS/DoA case, weighted Stansfield algorithm is shown to outperform other 
DoA fusion algorithms, such as maximum likelihood algorithm solved by iteration methods 
and linear least square algorithm [|2T1 . The weights for each DoA are their estimated error 
variances, which are function of RSS. The weighted Stansfield location estimate is given by 
ip,st = (A^W- x A S t) 



^A^W^st, where 



smw lj 



smW Nj 



COS(0i 



COSltfjvJ 



Xi sm(0ij — yi cos(#i, 



xi sm(V N ) - yi cos{V Nj 



(30) 
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and the weighting matrix is defined as W = diagprf, . . . , a 2 N ], where o 2 n is obtained by replacing 
true Do A 6 n with estimated Do A 6 n in © or ©. 

B. Simulation Settings 

We summarize the basic parameter settings in our simulations as follows. The CRs are 
uniformly placed in a circle of radius R = 150m and protective region R = 5m, centered 
around the PU with location ip = [0, 0] T for simplicity and without loss of generality. The PU 
transmit power P T is 20dBm (lOOmW) which is the typical radiation power of IEEE 802.11b/g 
wireless LAN transmitters for 20MHz channels in ISM bands. The measurement noise power 
Pm is — 80dBm (lOpW) which is a moderate estimate of noise introduced in the DoA estimation 
process, compare to the — lOOdBm (O.lpW) thermal noise power of the 802.11 WLAN channel. 
The path-loss exponent and shadowing variance are 7 = 5 and a s = 6dB respectively. The 
power and channel settings result in an averaged and minimum received SNR of lOdB and 
— lOdB, respectively. Each CR is equipped with an ULA with N a = 2 antennas, and uses 
N s = 50 samples for each localization period. The array orientation error with regard to the 
incoming DoA is assumed to follow uniform distribution given by 6 n ~ U{— 7r/3, 7t/3). We use 
the parameter 7r/3 for the orientation error since the maximum orientation error is 7r/2 which 
gives infinite error variance as can be verified through © or ©, thus n/3 is a practical value. 
Each data point in the figures is obtained from averages of 1000 CR placements (with random 
array orientation) and 2000 realization of RSS/DoA measurements if applicable. In the following 
subsections, we use these settings unless stated otherwise. 

C. Impact of Number of CRs 

We first study the impact of node density on the localization accuracy. The results for the 
RSS-only CRB (PT2)) and WCL algorithm are presented in Figffl It is observed that adding more 
CRs provides steady performance improvement for both RSS-only CRB and WCL algorithm, 
and the WCL algorithm has a constant RMSE gap of about 10m for sufficiently large number 
of CRs, say 40, compared with the RSS-only CRB. The results for the derived joint CRB using 
optimal DoA estimators (TT6l) and the weighted Stansfield algorithm are presented in Figfll The 
localization accuracy of the joint CRB and algorithm vary from 0.1 meters to 5 meters, compared 
with the accuracy range of 5 to 35 meters for RSS-only CRB and WCL in FigEfl indicating that 



13 



using both RSS and DoA measurements provides much better accuracy than using only RSS. 
The predicted accuracy by the joint CRB is achievable by the weighted Stansfiled algorithm, for 
large number of CRs. 

D. Impact of Channel Conditions 

We study the impact of channel conditions including shadowing variance and correlation 
distance on the localization accuracy. The considered values of shadowing variance vary from 4 
to lOdB, and the typical correlation distance for mobile communication is shown to be within 30 
meters ll22ll . The results for the RSS-only CRB and the WCL algorithms are shown in FigfJl It is 
observed that increasing shadowing standard derivation causes a linear increase in the RMSE, in 
that a 2dB increase of shadowing standard derivation results in 4 meters accuracy loss for both 
RSS-only CRB and the WCL algorithm. The performance loss caused by correlated shadowing 
is more significant for the RSS-only CRB than WCL. For example, when the correlation distance 
increases from meter to 30 meters under lOdB shadowing, the RMSE rises 5 meters for the 
RSS-only CRB, however the loss is about 2 meters for WCL. Therefore, WCL is robust to 
correlated shadowing. 

The impact of shadowing standard derivation on the accuracy of derived joint CRB, the 
joint CRB of lfT6ll and the weighted Stansfield algorithm are presented in FiglU Both results 
using optimal DoA estimator and the MUSIC algorithm are shown. We obtain a very surprising 
observation that when the shadowing standard derivation increases, the RMSEs of all considered 
cases decrease. This is because the RSS ip n is a log-normal random variable, when the shadowing 
standard derivation increases, the bell-shape PDF of the log-normal variable will become more 
spread-out, resulting in higher probabilities for larger RSS and lower RSS observations. Since 
the DoA estimation error variances © and © are inversely proportional to RSS, increasing 
shadowing results in more CRs with both good and bad DoA measurements. Therefore with 
proper weighting based on RSS information in the DoA fusion process, we may obtain better 
accuracy for large shadowing cases. The weighted Stansfield algorithm cannot achieve the joint 
CRB since it uses the estimated error variance for the weighting scheme; compared to the 
perfect knowledge of error variance obtained by the joint CRB. We can also observe that using 
MUSIC algorithm causes a 1 meter RMSE loss for weighted Stansfield algorithm, which is 
more significant compared to the joint CRB. Using MUSIC algorithm has small impact on 
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the joint CRB, since the increase in DoA estimation error is mitigated by making use of the 
perfect knowledge of each DoA's quality. Comparing the derived joint CRB and the CRB of 
|fT6l we conclude that the new CRB further reveals the potential achievable accuracy of the joint 
localization algorithms, therefore is more informative. 

The impact of correlation distance on the accuracy of derived joint CRB and the weighted 
Stansfield algorithm are presented in FigJ5J We observe that both the joint CRB and the weighted 
Stansfield algorithm are robust to correlated shadowing, in that the performance for joint CRB is 
almost constant for all X c values, and the RMSE loss for weighted Stansfiled algorithm is only 
0.2 meter from i.i.d shadowing to X c = 30m. The robustness of the joint CRB and weighted 
Stansfield algorithm to correlated shadowing is due to the fact that correlation only affect the 
RSS measurements which determine the variance of DoA measurements; however the Do As at 
each CR are corrupted by independent noise. Therefore the correlated shadowing does not have 
significant impact on the DoA fusion, and the resulting RMSE. 

E. Impact of Array Parameters 

In this subsection we study the impact of various array parameters, including number of 
antennas, number of samples and array orientation error, on the localization performance. The 
comparison results for varying number of antennas are shown in Fig|6] It is observed that the 
benefit of adding number of antennas from 2 to 3 is about twice compared with 3 to 4; and 
the benefit continues to saturate for adding even more antennas. This indicates that in practical 
systems, using a small number of antennas at each node is good from both hardware cost and 
localization accuracy point of view. 

The number of samples used in practical systems is affected by several factors such as signal 
bandwidth, sampling period, spectrum sensing regulations etc, which are beyond the scope of 
this paper. Therefore we provide results for some typical values of N s to gain some insights on 
the impact of this parameter on the RMSE, with results shown in FigjT] It is observed that the 
benefit of adding samples also tends to saturate when the number of samples is large enough, 
say 100. Comparing the results with Figj6] we conclude that increasing number of antennas 
introduces more significant improvement on RMSE. For example, a combination of 3 antennas 
and 50 samples achieves RMSE of 0.75 meters; however this accuracy is not achievable for 2 
antennas, even through we reach 150 samples. 
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The impact of array orientation error parameter T is studied by results shown in Figj8j It is 
observed that the joint CRB is robust to 8 T as the RMSE increases only 0.3 meters from 8 T = 
to 8 T = 5n/12. However the weighted Stansfield algorithm can only tolerate small orientations 
errors, and the RMSE rises by 1 meter from 8 T = tx/A to 8t = hix/12 due to difficulties to 
accurately estimate DoA estimation error variances for large orientation errors. 

F. Accuracy of Asymptotic CRBs 

We now evaluate the accuracy of the asymptotic CRBs for uniform random placement derived 
in Section IV compared with the exact CRBs conditioned on a specific placement derived 
in Section III. From (|22|) and (1281) we observe that the required number of CRs to make 
the asymptotic CRB accurate is a function of various array and channel parameters. It is 
straightforward to numerically investigate impact of each system parameter on the required 
N. Among all involved parameters, we notice that the value of Ro/R plays an important role in 
determining the required N. Thus we provide the comparison of exact CRBs after numerically 
averaging over uniform CR placements with the asymptotic CRBs for different R , with results 
shown in Figj9] and FigflO] for RSS-only bounds and joint bounds, respectively. For RSS-only 
bounds, the asymptotic CRB well approaches the exact bounds, in that the difference is less than 
1 meter for N = 40 when Ro/R = 0.1, and the accuracy of the asymptotic bounds is almost 
the same as the exact bounds for Ro/R = 0.3 and Ro/R = 0.5. Note that the RMSE increases 
as we increase Rq since we are excluding the nodes that are closer to the PU. On the contrary, 
since the variance of the asymptotic joint bound is larger than the asymptotic RSS-only bound, 
the approximation is less accurate. As we can observe from the figure, the approximation is 
satisfactory for N > 40 when Ro/R = 0.3m, and N > 30 when Ro/R = 0.5m. In practical 
systems, the circle with radius R correspond to the 'forbidden region' around the PU where 
activities of secondary network will severely interfere the PU and be interfered by the PU, the 
value of Ro depends on primary network requirements. However, if we consider the area ratio 
given by R^/R 2 , then both the asymptotic RSS-only and joint CRB provide satisfying accuracy 
for area ratio greater than 9%, which is a very reasonable forbidden region setting. 

We also evaluated the required number of CRs for asymptotic RSS-only and joint CRBs, 
predicted by (l22l) and (T28T) respectively, for varying Ro/R. The results are shown in Fig. [TT] 
Since it is not possible to directly relate the amount of deviation <5 in Theorem Q] and [2] with 
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the gap between the asymptotic CRBs and the exact CRBs, the choices of 5 and 77 become 

design parameters. As we have explained in Sec. [TV] (l22l) and (1281) tend to overestimate the 

required N, therefore we cannot select very small 5 and rj. In Fig. \TT\ S = |E F7 for 

L F 

RSS-only case, 5 — 2 || — E [F] || for the joint case and r/ = 0.15, which is the probability that 
the amount of deviation of the ensemble average is comparable to the statistical mean is 15%. 
We use a larger constant for the joint case since the variance of F is greater than that of F^, 
which results in a greater allowance on the deviation. Comparing the results in Fig. [H] with that 
in Fig. [9] and Fig. [TOl we observe that the predicted N for RSS-only case is accurate, in that at 
predicted N, the asymptotic RSS-only CRB in Fig. |9] well-approaches the exact CRB; on the 
other hand, the predicted TV for the joint case is in general more than enough. For example, 
when the R /R = 0.3, the predicted N is 90; however, result in Fig. [TO] indicates using 60 nodes 
already provides satisfying approximation accuracy. 

VI. Conclusion 

In this paper we present a novel analytical framework to evaluate the achievable performance 
of joint RSS/DoA algorithms for localization of PU transmitters in CR networks, which explicitly 
considers the dependency of DoA measurement quality on instant RSS. We first derive the joint 
CRB for a fixed CR placement, for both optimal DoA measurement and the MUSIC algorithm, 
respectively. We then derive the asymptotic joint CRB for uniform CR placement and provide 
theoretical and numerical study of its approximation accuracy. Using the derived joint CRBs 
in conjunction with numerical simulations, the impact of various system, antenna array and 
channel variables on the joint CRB and several practical localization algorithms is investigated 
and quantified. Based on the above analysis, we provide guidelines in terms of number of nodes 
and antenna array settings required to achieve certain localization accuracy, which are useful for 
practical localization system design. It is observed that the achievable RMSE decreases for large 
shadowing standard derivation due to cooperative diversity, and is robust to correlated shadowing. 
With a small number of antennas, say 2 to 3, and sufficient number of samples, jointly using 
RSS and DoA measurements reduces the RMSE by 90% compared with using only RSS. 
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Appendix A 
Derivation of Ft- and F 
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In this appendix we provide detailed derivations to obtain elements of F^ and Fg^- For F^, 
starting from (flOl) . each element in F^ is given as 
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where elements for all vectors of derivatives are given by 
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Applying the expectation in (I31H33I) with the derivatives results in (flTT) . 

For Fg^, starting from ([14)) . the second order derivatives of g n over xp is given by 
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Note that here we use the general form f((f) n ) to represent the function of RSS in the DoA 
estimation error variance, instead of using specific functions such as fcRB{4>n) or fMu(4>n)- 
Therefore the results in this appendix can be applied to derivation of Fg,^ of both assuming 
optimal DoA estimation or using MUSIC algorithm. When we take expectations of (|34l ), we 
apply the following property of taking expectation of a function f(x, y) of two random variables 
X and Y with known joint PDF p(x, y): E x ,y [f(%, y)] = Ey {Ejqy [f(x, y)]}, and get 
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Applying similar procedure to 
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Now we consider different /(0 n ) given by CRB or MUSIC, which results in 



E- 



f(<f>n) 



coP T e°s/M Pm 
dl N a 



p2 

— % 

/V 2 <t> 

a 



CRB, 
MUSIC, 



(37a) 
(37b) 



where the expectation in (I37bl) is an expectation of inverse of shifted log-normal variable, which 
is not directly obtainable and need to be calculated through numerical methods. 
The second order derivatives of h n over the PU coordinates are given by 
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where the second equality holds because the second order derivatives of h n are not related to 
random variables or <fi. Summation of corresponding elements in (|35l) . (l36l) and (1381) with 
some simplification will give us the elements of sub-FIM Fg^. 

Appendix B 
Proof of Theorem [Hand [2] 

We start by introducing the Chebyshev's inequality for random matrix, which will be the basic 
tool for the following proof. 

Theorem 3: (Chebyshev's Inequality for Random Matrix) Suppose X is a random matrix with 
mean E [X] and its second order moments exist. Then for any 5 > 0, we have 

E[||X-E[X]|| 2 ] 



Pr{||X-E[X]|| >5}< 
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where ||-|| denotes any matrix norm. 

Proof: From the Markov's inequality which states that if X is a nonnegative random variable, 
and E [X] is its mean, then for any 5 > 0, we have Pr{X > 5} < ^p-. Therefore, for scalar 



E[||X-E[X]|| 2 ] . 

-± — ji L . Taking square 



random variable ||X - E [X] f, we have Pr { ||X - E [X] f > 5 2 } < 
root of the expression on the left-hand-side finishes the proof. ■ 
To prove Theorem [Q we observe from the Chebyshev's inequality that we need to compute 
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As a result, we need to evaluate the following expression 
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The covariance of F-? appeared in the right-hand- side of (|40l) can be derived from ( TT9b as 
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From (|40k and (SB we obtain that 
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The application of Theorem [3] on 
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with result of (|42j) proves Theorem [IJ 



To prove Theorem [2l we observe from the Chebyshev's inequality that we need to compute 
E [||X - E [X] \\ 2 F ] for X = jfF, which requires evaluation of the following expression 
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F^ + Fg^ appeared on the right-hand-side of (|43l) can be derived as 
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where f n , E [f n ] and E [f 2 ] are given by (f27l) . We skip the derivations to obtain (l44l) which 
include tedious manipulations of ( TT9b and (|24|) . From (|43T) and (|44|) we obtain that 
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The application of Theorem [3] on II -^F — j? E [F] II with result of (1451 ) proves Theorem [2] 
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Fig. 1. Results for RMSE of RSS-only CRB and WCL algorithm with varying number of CRs, with uniform random placement 
in a circle with R — 150m. 
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Fig. 2. Results for RMSE of the joint CRB and weighted Stansfield algorithm with varying number of CRs, with uniform 
random placement in a circle with R = 150m. 
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Fig. 3. Results for RMSE of RSS-only CRB and WCL algorithm with varying shadowing standard derivation and correlation 
distance, with uniform random placement of 15 CRs. 




Fig. 4. Results for RMSE of the joint CRB, CRB of [16] and weighted Stansfield algorithm with varying shadowing standard 
derivation, with uniform random placement 15 CRs. 
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Fig. 5. Results for RMSE of the joint CRB and weighted Stansfield algorithm with varying correlation distance, with uniform 
random placement 15 CRs. 
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Fig. 6. Results for RMSE of the joint CRB and weighted Stansfield algorithm with varying number of antennas, with uniform 
random placement 15 CRs. For this scenario, N„ = 50. 
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Fig. 7. Results for RMSE of the joint CRB and weighted Stansfield algorithm with varying number of samples, with uniform 
random placement 15 CRs. For this scenario, N a = 2. 
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Fig. 8. Results for RMSE of the joint CRB and weighted Stansfield algorithm with varying array orientation error parameter, 
with uniform random placement 15 CRs. 
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Fig. 9. Results for exact and asymptotic RMSE of RSS-only CRB with varying number of CRs, for different Ro- 
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Fig. 10. Results for exact and asymptotic RMSE of the joint CRB with varying number of CRs, for different Rq. 
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TABLE I 

Summary of Variables Used in This Paper 



Symbol 


Meaning 


N, n 


number of CRs, index of CRs 


R, Ro 


radius of: PU coverage area, guard region 


p r iT 


2-dimensional coordinate of the n CR 


/> r iT 
l P = [ip.jp] 


2-dimensional coordinate of the PU 


£ P = [ip,yp] 


Estimated 2-dimensional coordinate of the PU 




distance from the n th CR to the PU, distance on x-axis, distance on y-axis 


Pt, Pa/, Pn 


PU transmit power, measurement noise power, received SNR 


co, 7 


channel multiplicative gain at reference distance, path-loss exponent 


S n , Xc 


shadowing variable, its variance, correlation distance 




Covariance matrix of: the PU location estimates, shadowing variables 




RSS measurement: in Watt, in dBm, mean in dBm 


d n , 9n 


incoming DoA and its estimate 


dn, d n , Ot 


ULA orientation, orientation error and its distribution parameter 


N s , N a , k 


ULA parameters: number of samples, number of antennas, array constant 


2 2 
a n,CRB' a n,MU 


DoA error variance of the n th CR given by: CRB, MUSIC algorithm 


F 5,> F §is>' F 


FIM for: RSS-only, DoA given RSS, joint RSS/DoA scenarios 


a, P, e 


supporting constants: a 4 aa^g/g^ p = ^ffi^-i)' 6 = 


So, V 


deviation distance and corresponding probability used in Theorem Q] and [2] 


s, cj>, 


vector form of: shadowing variables, RSS in dBm, DoA estimates 



